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We have probed the Landau levels of nodal quasi-particles by tunneling along a nodal direction 
of (110) oriented YBa2Cu307_ a; thin films with a magnetic field applied perpendicular to the CuOi 
planes, and parallel to the film's surface. In optimally doped films and at low temperature, finite 
energy nodal states are clearly observed in films thinner than the London penetration depth. Above 
a well defined temperature the order parameter reverts to a pure rf-wave symmetry. 
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The pairing interaction responsible for superconductiv- 
ity in the high temperature superconductors is still under 
debate. There is much experimental evidence showing 
that the order parameter has a dominant d x 2_ y 2 symme- 
try Q) ) but important questions such as the nature of 
quasi-particle states in the under-doped region and the 
existence of a minority component of the order parameter 
are still under debate 0|. Studying the density of states 
when a magnetic field is applied along the c-axis can lend 
important information on both of these questions. 

Gor'kov and Schrieffer [J] have remarked that in a 
d-wave superconductor nodal quasi-particles undergo in 
their orbital motion a series of Andreev - Saint-James 
reflections from nearby lobes of the order parameter, re- 
sulting in finite energy states. This states, as noted by 
Anderson [Bj], carry a current around the Fermi surface. 
Each time an electron, say, is reflected as a hole, the 
missing pair reappears as a Cooper pair in the conden- 
sate. Such currents can be described by an id xy com- 
ponent which, as shown by Laughlin [f|, lowers the free 
energy of the superconducting state under the applied 
field because of the moment they produce. Minimization 
of the free energy with respect to the amplitude, 5 xy , of 
the minority component, lends at T = the law: S xy — 

Here v = y/viv?, where v\ is 



aH 1 / 2 with a = hvJdz. 

the Fermi velocity, vi = j-^- at the node direction, and 
A is the main d x 2_ y 2 superconducting gap. A weak first 
order phase transition to a pure d-wave symmetry was 
predicted to occur at a temperature UbTcf = b5 xy (H), 
where b is a universal constant equal to 0.52. It was 
however pointed out that in the mixed state Meissner 
currents would cause a large Doppler shift of the nodal 
states energy, making their observation difficult if not 
impossible [7j. 

In this letter we show how this difficulty can be re- 
solved by taking measurements on films thinner than the 
London penetration depth, and conclude that a field in- 
duced minority component exists below a well defined 
temperature Tcf- Our results are in semi-quantitative 
agreement with Laughlin's theory [f| with proportional- 



ity constant, b = 0.34, smaller than the predicted one. 

One of the most powerful tools for probing the order- 
parameter is tunneling spectroscopy. For conventional 
superconductors, the tunneling spectrum is directly pro- 
portional to the superconductor's density of state and 
exhibits two peaks at biases corresponding to ±A, where 
A is the superconducting gap [![. For a c£-wave order- 
parameter symmetry, zero energy surface states are 
formed when quasi-particles undergo successive Andreev 
- Saint-James reflections from order parameters whose 
phase differ by tt, as occurs upon reflection at a surface 
perpendicular to a nodal direction. Tunneling along a 
node direction, i.e. the (110) direction, shows a zero bias 
peak in the tunneling spectrum (9J due to these states 
[lfj . An additional minority id xy component would show 
up as peaks of the junction's conductance at biases equal 
to S xy , replacing the zero bias peak [§]. 

According to theoretical predictions 0, @ , tunnel- 
ing spectral peaks should thus appear when the field is 
applied perpendicular to the Cu02 planes in films hav- 
ing a nodal orientation. The presence of such peaks was 
not ed p reviously, but received a different interpretation 
[ll , TJ H| , which is that Meissner screening currents will 
Doppler shift the zero energy surface states. It was shown 
[14j that the split of the zero bias conductance peak is al- 
ways larger in increasing fields than in decreasing ones, a 
behavior that can be understood on the basis of the hys- 
teretic behavior of screening currents. These currents are 
of two kinds: Meissner currents on the scale of the Lon- 
don penetration depth, and Bean currents due to vortex 
pinning. Meissner currents are large in increasing fields, 
due to the Bean Livingston barrier that retards the pene- 
tration of vortices, and small in decreasing fields because 
there is no barrier against vortex exit. As for Bean cur- 
rents, they reverse sign upon field reversal, and are typ- 
ically weak compared to Meissner currents in increasing 
fields. In experiments performed on YBa2Cu307_ a: thick 
films having the (110) orientation, Beck et al. [IH have 
reported on finite bias conductance peaks that follow the 
H 1 / 2 law in decreasing fields where screening currents are 
weak, and attributed them to a field induced id xy compo- 
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FIG. 1: A scheme of the samples' configuration. The samples' 
surface is in the (110) direction, and one of the surface edges 
is parallel to the c-axis. The magnetic field is applied parallel 
to the surface and to the c-axis. 



nent. This conclusion has however remained controver- 
sial as the intensity of screening currents is not known 
with certainty. 

In an effort to establish more firmly the origin of the 
field induced finite bias conductance peaks, we have made 
a series of tunneling measurements as a function of field 
and temperature on films thinner than the London pen- 
etration depth. Both the Bean Livingston and the Bean 
currents must then be quite weak, since screening cur- 
rents generated on opposite faces of the film cancel each 
other. If the origin of the finite bias peaks seen in thick 
films is a Doppler shift of zero energy surface bound 
states, they should not occur at all in thin films, either in 
increasing or in decreasing fields. If, on the other hand, 
they are basically due to field induced finite energy nodal 
states, they should persist in the thin films with no field 
hysteresis, and follow the if 1 / 2 . As we now describe in 
detail, this is indeed what we have found. 

We fabricated (110) oriented YBa2Cu307_ x thin films 
by DC off-axis sputtering on (110) SiTrOs substrates 
with one of the surface edges parallel to the c-axis. A 
PrBa2Cu307_a; buffer layer was first deposited by off- 
axis RF sputtering in order to reduce the growth of (103) 
oriented grains [16j . The thickness of the samples is less 
than 500A. X -ray diffraction reveals peaks corresponding 
to the (110) orientation, the Tc of the samples ranged 
from 84 K to 89 K down-set. The films' surface is rela- 
tively smooth, having an average surface roughness of few 
tens of angstroms according to atomic force microscope 
measurements. We verified the in-plane orientation by 
comparing the normal state resistance along the two dif- 
ferent directions and found the expected anisotropy [l6| . 

The tunneling junctions are produced by pressing In- 
dium thin pads on fresh samples, as described elsewhere 
14 , Il7| . The configuration of the sample and the field 
orientation is shown in Fig. [T] The Indium-Oxide layer, 
which is created in the area of contact, results in junc- 
tion resistances in the range of 1-15 fl. This process 
may reduce the oxygen concentration in the junctions' 
area, especially in very thin samples. In order to alle- 
viate this effect, we grew our films in the presence of 
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FIG. 2: (color online) (a) The conductance spectrum of a thin 
sample obtained when increasing and decreasing the magnetic 
field. Full black line: H=0 T. Arrows facing up(down) repre- 
sent increasing(decreasing) field, (b) The conductance spec- 
trum obtained on a thick sample. Full black line: H=0, thin 
dashed line: H=4 T increasing field, thick dashed line: H=4T 
decreasing field, (c) Spectral peak value versus the square root 
of applied field in decreasing fields measured on a thin sam- 
ple at 4.2 K. The solid line is a linear fit with a slope of 
0.65 mV/T 1 ' 2 . We note that the slope found here is slightly 
smaller than that found for thicker films (see text). 



a reduced vapor pressure of water, which results in in- 
creased oxygen content in the film [l8j |. The junctions 
are stable under thermal cycles, and the junctions' re- 
sistance remains unchanged if kept in a helium gas en- 
vironment. This allows us to measure the temperature 
dependence of their tunneling characteristic, if necessary 
by performing successive cooling runs under various ap- 
plied fields. The planar junction technique averages the 
tunneling conductance at a macroscopic length scale in 
contrast to scanning tunneling microscopy measurements 
which give sharper peak features [19(. However, its sta- 
bility over thermal and magnetic fields cycles is of great 
advantage for the present work. 

In Fig. [2k we show the differential conductance ob- 
tained on a thin film at 4.2 K, measured at magnetic fields 
H—Q, 6 and 12 T, in increasing and decreasing fields. At 
zero magnetic field, there is a clear peak at zero bias, as 
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expected for a (110) oriented film [9(. This peak splits 
into two spectral peaks in the presence of a finite mag- 
netic field. The conductance spectra in increasing and 
decreasing fields are almost identical as opposed to the 
high hysteretic behavior of the spectral peaks in thicker 
samples as shown in Fig. (2b and in Ref. FJ, ll7[ . This 
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highly hysteretic behavior has been ascribed to strong 
screening currents in increasing fields and much weaker 
ones in decreasing fields as discussed above. IU 13| ■ The 
negligible hysteresis observed in our thin samples implies 
that the screening currents are very low even in increas- 
ing fields and that the two finite bias peaks cannot be due 
to a Doppler shift of the zero energy bound states. The 
difference between the behavior of thin and thick films 
is particularly dramatic in increasing fields, where the 
spectral peaks cannot be identified anymore at all in thick 
samples. This is because the Doppler shift, due to Meiss- 
ner currents, widens the peaks and pushes them into high 
bias where they are merged with the background and the 
main gap structure (Fig. Wp)- 

We define the minority gap peak value, 8, as half the 
distance between the positive and negative bias conduc- 
tance peaks (see Fig. [2p). In our thin films, S(H) follows 
the H 1 / 2 behavior as shown in Fig. [2fc. This is in agree- 
ment with measurements performed earlier on thick films 
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in decreasing fields 



20( and confirms the assumption 
that Meissner currents are weak in that case, but may be 
not negligible in view of the smaller slope obtained here. 

Measuring the tunneling spectrum at various magnetic 
fields and temperatures, we find that in a fixed magnetic 
field, the minority gap peak value disappears above a well 
defined temperature. For example, we show in Fig. [3^ 
the measured conductance in a field of 7 T at various 
temperatures. At 1.6 K, there are two clear peaks at ± 
1.75 mV. As the temperature rises, the conductance at 
zero bias increases until the finite bias peaks completely 
disappear at 8.5 K. The conductance at high biases is 
independent of temperature as expected for tunneling, 
assuring us that the junction characteristics remain in- 
tact. 

However, in order to be able to properly assert if there 
is a temperature induced modification of the density of 
states at low bias, one must take into account the effect 
of thermal smearing on the tunneling spectra 0] . There- 
fore, we have convoluted the tunneling spectra measured 
at the lowest temperature (below 1.6 K) with the deriva- 
tive of the Fermi-Dirac distribution function at the de- 
sired temperature. Previous scanning tunneling measure- 
ments showed that the density of states hardly changes 
at such low temperatures [2l|. Therefore, using our pro- 
cedure, we have calculated a thermally smeared curve 
for higher temperatures and thus estimated the effect of 
thermal smearing on the low energy spectrum. Figure [3b 
shows the measurement taken at 8.5 K (red full line) and 
the result of smearing the 1.6 K measurement to 8.5 K 
(black dashed line). While the smeared curve still shows 
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FIG. 3: (color online) (a) Conductance measured under a field 
of 7 T at various temperatures, (b) Comparison between the 
measured data at T=8.5 K, H—7 T(full line) and the result 
of smearing the data measured at T=1.6 K, H—7 T to T=8.5 
K (dashed line). 



two distinct peaks, the measurement does not, but rather 
exhibits a peak at zero bias. It is therefore clear that 
temperature induces a change in the density of states. 
At biases higher than 5 mV, the two curves are identical. 

A representative graph of the minority gap value ver- 
sus temperature for various values of the applied field 
is shown in Fig. The minority gap value is ap- 

proximately constant at low temperatures and disappears 
abruptly at a temperature, which we define as Tcf- The 
minor initial rise in the sub gap peaks' value is an effect 
of thermal smearing. 

The value of Tqf increases with the applied field 
and the resulting low temperature conductance peak 
value. We find that they obey the universal linear re- 
lation ksTcF — 0.345 as shown in Fig. [4)d. The data 
points in the graph were taken from five different sam- 
ples with different characteristics (one of the points was 
measured on a thick sample in decreasing field), which 
implies the universality of the proportionality constant. 
The measured slope is smaller than Laughlins' prediction 
^bTcf = 0.525. Laughlins' calculation was performed 
for bulk materials, while our measurements were done 
on very thin films, which are highly susceptible to sur- 
face effects. Laughlin did not mention the length scale of 
the currents responsible for the dipole moments. If the 
thickness of the film is comparable to the aforementioned 
length scale, interaction between the currents resulting in 
a partial cancelation of the opposite currents should be 
considered. This might contribute to the disagreement 
between the theoretical ratio and our experimental re- 
sult. 

A way to distinguish between the Doppler shift effect 
and the minority component induced by the applied field, 
is by their respective temperature dependence. A first 
order phase transition that leads to the abrupt disap- 
pearance of the minority orderparameter at high tem- 
peratures has been predicted |6|. On the other hand, 
screening currents and the accompanying Doppler shift 
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FIG. 4: (color online) (a) The spectral peak values versus 
temperature obtained from different applied magnetic fields, 
(b) The temperature at which the spectral peak disappears 
versus the peak value at low temperatures. 



superfiuid currents are strong (and at the origin of the 
square root dependence) and predicted to render the mi- 
nority component unobservable 0]. 
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should persist at high temperatures since the superfiuid 
momentum is independent of temperature at tempera- 
tures significantly lower than Tq [221 ] . In addition to the 
small field hysteresis in the position of the finite bias 
peaks, their sharp disappearance at a well defined tem- 
perature rules out that a Doppler shift of the zero energy 
surface bound states is at their origin. 

The proposed interpretation of the experiments that 
we have described requires that the scattering time of 
nodal quasi-particles be long compared to the time it 
takes them to complete a Saint- James cycle. This ap- 
pears to be true at optimum doping, but not away from 
it as we show in a separate publication (23j | . 

In conclusion, we claim that the field induced spec- 
tral peaks seen in nodal tunneling in very thin films re- 
sult from the formation of a minority order parameter, 
and cannot be attributed to a Doppler shift of the zero 
energy states due to Meissner screening currents. The 
square root dependency of the spectral peak positions in 
the magnetic field strength, the absence of field hysteresis 
and their abrupt disappearance with temperature are all 
in favor a field induced id xy minority imaginary order pa- 
rameter at low temperatures and a transition to a pure d- 
wave order parameter at high temperatures, as predicted 
by Laughlin. Our findings are in general agreement with 
his model, but a comparison of the way in which the tun- 
neling density of states evolves with temperature would 
require a more detailed theory than is now available. The 
absence of a minority order parameter inferred from the 
heat capacity square root field dependence [24| is not in 
contradiction with our conclusions, because these mea- 
surements were performed on much thicker crystals where 
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